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ABSTRACT 1 

A numerical analysis was perCormed to compare natural convection 
velocities in two-dimensional enclosures of various shape. The following 
shapes were investigated; circle, square, horizontal and upright 2x1 | 

aspect ratio rectangles, horizontal and upright half-circles, diamond | 

(square oriented with diagonal vertical) and triangle (equilateral with \ 

horizontal base). In all cases, the length scale in the various dimen- | 

sionless parameters, such as Rayleigh number, is defined as the diameter of ! 

it 

equal area circle. Natural convection velocities were calculated lur I 

Rayleigh numbers of 1000 and 5000 with the temperature difference taken to | 

be across (a) the maximum horizontal dimension, (b) the median horizontal 
line (line through centroid) and (c) the horizontal distance such that tltc 
temperature gradient is the same for shapes of equal area. A Rayleigh s 

number of 1000 is within the "low Rayleigh number" range: for agreement with 
first order theory for circular enclosures. A Rayleigh number of 5000 is 
slightly out of this range. For the class of shapes including the square, 
upright half-circle and upright rectangle, the computed velocities were | 

found to agree very closely with that of the equal area circle when the I 

temperature difference is taken to be across the maximum horizontal {| 

dimension (condition (a)). The velocities for the horizontal rectangle and ll 

half circle were found to be approximately one-half that of the equal area | 

circle for the same condition. Better overall agreement among all shapes 
was obtained by setting the temperature difference across a distance such 
that the temperature gradients were equal for shapes of equal area. 
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NOMENCLATURE 


Derr^ription 

equal area circle diameter 
gravity force 
circle radius 

3 

Rayleigh number = — 

temperature difference 
time 


velocity 

thermal diffusivity 

volumetric coefficient of thermal expansion 

dynamic viscosity 

kinematic viscosity, y /p 

stream function 

density 
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1. INTRODUCTION 

A number of fluid mechanics experiments are to be performed aboard 
orbiting spacecraft, primarily to investigate the possible use of the 
near-zero gravity environment in various materials processing applications. 

A knowledge of the intensity of convective stirring due to residual accele- 
rations is required in order to properly plan and evaluate the experiments. 
Various shaped containers will be used in these experiments. Theoretical 
results have been obtained which yield exact predictions ol natural convec- 
tion velocities in an idealized container, i>e., the two-dimensional cir- 
cular enclosure (Ref. 1), It would be helpful for estimation purposes if a 
reasonable means existed for the extrapolation of the circular enclosure 
results to more complex shapes. A previous study demonstrated extremely 
good agreement in computed natural convection velocities for circular acid 
square enclosures of equal cross-sectional area (Ref. 2). Dressier (Ref. !3) 
made use of this noted agreement in circular and square enclosure results, 
in addition to the more general results documented herein, to analyze 
natural convection in the proposed Lai, Kroes and Wilcox crystal growth 
experiment to be performed Infight on Spacelab 3. The purpose of this 
investigation is to develop reasonable extrapolation criteria by comparing 
numerically computed natural convection velocities for various two- 
dimensional enclosure shapes with the circular enclosure results. The 
Lockheed-developed General Interpolant Method (GIM) computer code (Ref. 4) 
was used in the numerical computations. 
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2. METHOD OF APPROACH AND NUMERICAL SIMULATION 

The following set of two-dimensional enclosure shapes was considered; 
circle, square, 2x1 aspect ratio rectangle in both upright and horizontal 
orientations, half-circle in both upright and horizontal orientations, 
diamond (square oriented with diagonal vertical) and triangle (equilateral 
with horizontal base). The baseline condition was selected for equal area 
with the temperature difference set across the maximum horizontal dimen- 
sion. The initial temperature distribution was based on a uniform hori- 
zontal gradient with the boundary points held constant in time. The gravity 
vector was considered to be constant in the downward direction. The fluid 
was assumed to behave as a boussinesq fluid in its thermal expansion 
characteristics . 

The numerical simulation was based on dividing the various enclosures 
into a computational grid consisting of a network of generalized quadri- 
lateral elements with curvilinear sides, with the nodal points located at 
the four corners of each element. Each of the enclosures considered in this 
study was treated as a generalized quadrilateral region divided by inter- 
polation into an array of 20 x 20 elements. For example, the circle was 
treated as a four-sided figure, each side being a quarter-circle arc. 

The triangle was treated by forming a parallelogram such that two 
adjacent sides and the diagonal formed the desired triangle. The entire 
parallelogram was divided into a computational grid, but the diagonal points 
were treated as boundary points. The parallelogram with the diagonal 
boundary thus formed two independent triangular regions. The computational 
grids for the enclosures considered in this study are shown in Figs. 1-5. 
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The grids for the square and diamond are the same since the diamond is 
simply a square with the diagonal in the upright position. The rectangle 
and half-circle grids are identical in either the upright or horizontal 
positions* 
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3. RESULTS 

The velocity histories for the various shaped enclosures are shown in 
Figs. 6 and 7 for I^yleigh numbers of 1000 and 5000. The various enclosures 
are all of equal area, with the diroensxonless velocities and Rayleigh 
numbers defined based on a temperature difference ^dT across the maximum 
horizontal dimension, and a length scale equal to the diameter d of the 
equal area circle. Note that, for a Rayleigh number of 1000 the steady 
state velocities of the equal area circle and square are very close to 
identical, with the upright rectangle and half-circle being within about 5 % 
of the circle value. The horizontal rectangle and half-circle, and the 
triangle steady state velocities form another grouping of values approx- 
imately one-half the circle value. The diamond results appear approximately 
midway between the two extremes. Roughly the same trend of steady state 
values appear for a Rayleigh number of 5000. Overall, the steady state 
velocities for the various equal area enclosures are in agreement within a 
factor of approximately 2. 

Response times for the various enclosure results are shown in Fig. 8 
for both 1000 and 5000 Rayleigh number. Agreement within a factor of 
approximately 2 is shown for all enclosures. These response times are 
defined as the time required to reach the fraction 1 - 1/e (0.632) of the 
maximum or steady state velocities, whichever is greater, in Figs. 6 and 7. 

An attempt was made to find other bases for correlating the data to 
yield better agreement between the results for the various shaped enclo- 
sures. Including the base line correlation, described earlier, correlations 
were made for the following seta of conditions concerning the distance over 
which the temperature difference was taken: 

o Temperature difference across maximum horizontal dimension 
(base line correlation). 
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0 Tamperatrure difference across median horizontal line (line 
through centroid). Only the triangle and horizontal 
half "Circle results changed in this correlation from the 
baseline correlation. 

0 Temperature difference taken across a distance such that the 
temperature gradient is the same for all shapes. 

The steady state results for the above sets of correlations are sum- 
marized in Fig. 9 for a Rayleigh number of 1000, and in Fig. 10 for a 
Rayleigh number of 5000. For each of the correlation sets, the length scale 
in the Rayleigh number and dimensionless velocity is taken to be the diam- 
eter of the equal area circle. A comparison of the three corrulation sets 
shows improved correlation for the triangle and horizontal half-circle by 
taking the temperature dif Terence across the median horizontal line rather 
than the maximum horizontal dimension. Generally better correlation is 
obtained by setting the cemperature gradients equal. This is equivalent to 
taking the temperature difference across a distance equal to the diameter of 
the equal area circle. 

Computer generated streamline, absolute velocity and temperature con- 
tour plots for all of the enclosures are shown in Figs, 11 through 34 at 
steady state for both 1000 and 5000 Rayleigh numbers. 
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4, CONCLUSIONS 

Natural convection velocities withxn two~dimenoional enclosures of 
various shape, for Rayleigh numbers up to at least 5000, may be estimated 
with reasonable accuracy by considering the area to be equivalent to a 
circle of equal area. For the class of figures including the square, 
upright half-circle, and upright 2 x 1 aspect ratio rectangle, excellent 
agreement is obtained by considering the' temperature difference across the 
maximum horizontal dimension to be equal to that across the equal area 
circle horizontal diameter, the agreement of the two upright oblong shapes 
indicates probable agreement for any similar upright oblong shape of roughly 
the same aspect ratio, 

Tlie horizontal oblong shapes have natural convection velocities 
approximately one-half that estimated based on the equal area circle with 
the temperature difference across the maximum horizontal dimension taken to 
be equal to that across the equal area circle diameter. As with the upright 
shapes, this probably indicates a general pattern for similar oblong shapes 
of roughly the same aspect ratio. 

Better overall agreement among all shapes is obtained by setting the 
temperature difference across a distance equal to the diameter of the equal 
area circle, thus making the temperature gradients equal for any shape. 
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Fig. 5 - Computational Grid: Triangle 
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Velocity Histories for Natural Convection in Two-Didensional Enclosures 
of Variotis Shape and Rayleigh Number of 5000 
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Fig. 10 - Comparison of Steady State Velocities for Natural Convection 
in Two-Dimensional Enclosures of Various Shape and Rayleigh 
Number of 5000 






Absolute Velocity Contour Plots for Circular Enclosure 





Fig. 13 - Temperature Contour Plots for Circular Enclosure 




Streamline Plots for Square Enclosure 
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Absolute Velocity Conx.,ur Plots for Square Enclosure 
















Ratio Rectangular Enclosure 









Streamline Plots 






Aspect Ratio 




Temperature Contour Plots for Horizontal 2x1 Aspect Ratio Rectangular Enclosure 
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Absolute Velocity Contour Plots for Diamond Enclosure 






Triattgular Enclosure 
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Eig. 33 — Absolute Velocity Contour Plots for Triangular Enclosure 



LMSC-HREC TR D867666 


OF POOR QUALITY 



39 

LOCKHEED-HUNTSVILLE RESEARCH & ENGINEERING CENTER 








Fig. 34 - Temperature Contour Plots for Triangular Enclosure 


